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We have performed extensive simulations of random sequential adsorption and diffusion of k- 
mers, up to k = 5 in two dimensions with particular attention to the case k = 2. We focus on 
the behaviour of the coverage and of vacancy dynamics as a function of time. We observe that 
for k = 2, 3 a complete coverage of the lattice is never reached, because of the existence of frozen 
configurations that prevent isolated vacancies in the lattice to join. From this result we argue that 
complete coverage is never attained for any value of k. The long time behaviour of the coverage is not 
mean field and non analytic, with f -1 / 2 as leading term. Long time coverage regimes are independent 
of the initial conditions while strongly depend on the diffusion probability and deposition rate and, 
in particular, different values of these parameters lead to different final values of the coverage. The 
geometrical complexity of these systems is also highlighted through an investigation of the vacancy 
population dynamics. 
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I. INTRODUCTION 

Random Sequential Adsorption (RSA) of hard parti- 
cles on lattice is employed to describe a wide class of 
irreversible phenomena in physics, chemistry and biol- 
ogyB The prototype oLthis process is the dimerization of 
monomers on a chain.EI Each monomer occupies initially 
a site of the chain and can react irreversibly with one of 
its neighbors to form a dimer: irreversible dimerization 
takes place sequentially until no more pairs of neighbor 
monomers are available. According to FloryS this pro- 
cess can be viewed as the random sequential adsorption 
of non-overlapping dimers onto a one dimensional lattice: 
a dimer can be irreversibly deposited if two neighbor free 
sites are available. At a given stage only single isolated 
sites will be available and further deposition will be no 
more allowed. At that stage, the density of dimers is fi- 
nite and takes a non random value in the limit of infinite 
system, the so called jamming coverage 9j. This study 
was followed by the introduction and investigation of a 
variety of RSA,-processes involving different objects, lat- 
tices and rulesE The quantities of interest are in general 
the time behavior of the coverage 9(t) and its jamming 
value. Exact solutions are available in few cases, typically 
for dimers in one dimension, whereas only numerical or 
approximate results .exist for more complex objects and 
in higher dimension, aau 

In recent years, the class of systems for which RSA 



yields effective description has widened by the observa- 
tion that in many experimental processes particles also 
diffuse on the lattice after deposition. Diffusional Relax- 
ation (DR) affects the spatial distribution of free sites, al- 
lowing for further deposition of particles whenever neigh- 
bor empty sites are created. This produces a different 
time behavior of 9(t) and different values of the asymp- 
totic coverage 9(oo). The effect of DR on iRSA models 
has been investigated by many authors (seaffl and Refs. 
therein), often displaying the emergency of non trivial 
behaviors both for 9(t) and 9(oo). In the case of dimers 
on a line, for instance, the introduction of DR produces 
a very slow increase of coverage with time, terminating 
eventually with the complete fillingj-^Ioo) = 1. Exact re- 
sultsfl □ and numerical simulationaid'Lil show that, after 
a transient depending on the initial conditions, 1 — 9(t) 
decreases asymptotically as i -1 / 2 , and thus an infinite 
time is needed to reach the completely filled state. The 
exponent 1/2 is due to the relevance of statistical fluc- 
tuations in the formation of neighbors free sites, so that 
any computation based on mean field treatments is not 
effective. It is also for this reason that solutions of RSA 
problems with Diffusion (RSAD) are very race and usu- 
ally limited to one dimensional systems (sedi3 and Refs. 
therein) and some particular casesE'u Some Iwe dimen- 
sional models have been recently investigatedoli3ii3, and 
features sometimes different from the one dimensional 
case have been observed. Eisenberg and Baram have 
found a stable asymptotic coverage 9(oo) < 1 for the 
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RSAD of crosses on a square latticeE! whereas-such a be- 
havior is not observed in the case of squarescj, in which 
the complete coverage is eventually reached for an infinite 
system, although through coexistence of different phases. 
The pure RSA of dimers in two. dimensions was first in- 
vestigated by Nord and EvansE3 and was latec-extended 
also considering anisotropic deposition rates .EJ Grigera 
et al.t3 computed the asymptotic coverage 6(00) in pres- 
ence of diffusion and rotation of dimers and found that 
it is less than unity. They also observed indications for 
an exponential approach to the limit coverage 9 (00). 



and of their dependence on the deposition and diffusion 
rates. In Sec. [V the results are put in relation with the 
behaviour of the vacancies in the lattice, whose statistical 
behaviour is also studied. The fc— mers RSAD model with 
2 < fc < 5 is addressed in Sec. |v|, where similarities and 
differences with the case k = 2 are pointed out. Finally 
sec. VI is devoted to a summary of our main findings and 



to some general conclusions. 



II. MODEL AND COMPUTATIONAL DETAILS 
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FIG. 1. Temporal behavior of density of empty sites for 
different sizes of the lattice: L = 40 (solid line), L — 100 
(dotted line), L = 120 (dashed line), L = 140 (long dashed 
line). The initial condition is 8(0) = 8j. 

In this paper we present extensive Monte Carlo (MC) 
simulations of the RSAD of linear particles on a square 
lattice. In the systems under consideration hard par- 
ticles, the A:— mers, are dropped on an initially empty 
square lattice and are allowed to diffuse. Each fc— mcr 
occupies k consecutive sites and cannot overlap to the 
others. We consider different initial conditions as well 
as a set of different deposition and diffusion rates and 
probabilities. Whereas the former only have influence on 
initial and intermediate stages of the coverage, the latter 
influence the late stage of dynamics and different values 
may lead to different final states. In fact, we show that 
the value of the final coverage depends on these parame- 
ter and, as expected, decreases with increasing fc. In this 
paper the focus will be mainly on the case fc = 2, since 
simulations show that its main features are also shared 
by systems with higher values of fc, while the time re- 
quired for simulations rapidly increases with fc. For this 
reason Sees. are dedicated mainly to the RSAD of 

dimers. 

In Sec. |lj we define the model and the procedure for the 
Montecarlo simulation. We also test possible algorithm 
to save computational time when simulating such slow 
relaxation processes, and describe the one we adopted. 
Section [II is devoted to investigate the dynamical fea- 



A. RSAD dynamics 

The process that we have more extensively investigated 
in this paper consists in the deposition and diffusion of 
dimers on a square lattice containing N — L x L sites 
with periodic boundary conditions. Each dimcr occupies 
two nearest neighbor sites, or equivalently a bond, and 
cannot overlap with other dimers. For t > it is free 
to diffuse along the direction of its axis towards empty 
sites. We have assumed deposition to occur with equal 
probability on horizontal and vertical bonds and we have 
considered two possible different initial conditions: a) 
jamming coverage, i.e. 8(0) = 8j and b) empty lattice, 
i.e. 8(0) = 0. In the case a) the initial stage consists of 
pure RSA. Dimers are deposited on the lattice until no 
more free pair of sites is available and the jamming state 
is reached. The jamming coverage for this kind of model 
is known from previous workliU to be 8j = 0.90654. It is 
reached exponentially and is higher than the correspond- 
ing one in one dimension. At this stage deposited dimers 
are allowed to diffuse and more dimers can be possibly 
added when suitable space is made available. In case 
b) the deposited dimers can move from the beginning, 
so that diffusion in the initial stages largely dominates 
the deposition process, giving rise to a large amount of 
vacancies. 

TABLE I. Coefficients of fit Eq. (f| for 8(0) = 9j. The 
curve for p = 0.5 can be adequately fitted only retaining the 
terms up to the third order. 
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ri 


?~3 
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0.2 
0.5 
0.9 


0.00150 
0.00110 
0.00120 


0.040 
0.104 
0.252 


2.167 
0.607 
0.702 


7.447 
-1.105 
-3.983 


7.456 
4.799 


TABLE II. Coefficients of fit Eq. (g) for 6(0) = 0. 
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r3 
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0.2 
0.5 
0.9 


0.0016 
0.00084 
-0.0014 


0.0628 
0.1239 
0.3797 


1.358 
0.4223 
-0.6506 


2.152 
0.2943 
1.021 


-4.836 
-0.3776 
-0.4146 



tures of the coverage process as well as of the final state 
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FIG. 2. Temporal behavior of density of empty sites (solid 
lines) for p = 0.2 (a), p = 0.5 (b) and p = 0.9 (c). The long 
dashed lines are fits to the data according to Eq. (^). The 
initial condition is 8(0) = 8j and the size of the system is 
L = 100. 

We have computed the time behavior of the coverage 
and of the abundance of different kinds of vacancies in 
the above processes by means of MC simulations. At 
each time step deposition of dimers is attempted with 
probability p and diffusion with probability 1 — p. Both 
directions of diffusion are equiprobable. Deposition hap- 
pens at rate R, i.e. pRN dimers are tried to be de- 
posited on the average in a MC step and a unitary jump- 
ing rate d = 1 follows for diffusion by assuming unitary 
lattice spacing and simulation (micro) time discretization 
5t = l/N. 



B. Saving time algorithms 



The above dynamics is extremely time consuming to 
simulate, mainly in the late stage of the process, since 
due to the exclusion rules most of diffusion and depo- 
sition attempts fail. A large amount of computational 
time can however be saved and dynamics sped up by 
suitable algorithms that maximize the number of suc- 
cessful events. We have adopted the strategy of storing 
in memory the dimers that can successfully diffuse along 
a given direction, ridifH and n d ifv and the bonds avail- 
able for deposition ridepH and n dep y respectively for the 
horizontal and vertical directions. 
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FIG. 3. Temporal behavior of density of empty sites (solid 
lines) for p = 0.2 (a), p = 0.5 (b) and p = 0.9 (c). The 
long dashed lines are fits to the data according to Eq. (^J). 
The initial condition is #(0) = and the size of the system is 
L = 100. 

By accounting for the different probabilities and rates 
at which these actions may take place, it turns out that 
on the average a successful event occurs with probability 

_ / n depH + n depV \ l—p f ridj-fH + n dl f V 
Psucc-py m )+ 2 \ 2N 

(1) 

Since the law of large numbers is expected to hold (at 
least for system size large enough), at each step of the 
simulation the MC time is increased by the amount 



At = 



1 



Ps 



(2) 



We adopted the above method that produces a time sav- 
ing of about one order of magnitude with respect to the 
usual MC simulation with equivalent results. 

As a further check we [have also tested the method 
adopted by Gan and WangtLS, which is based on the prob- 
ability distribution of time lags At between two consec- 
utive successful events 



P(At 



Ps 



^(l Psucc) 



m — 1 



(3) 



and also yields results equivalent to those of usual MC. 
The identity of results confirms the reliability of both of 



these methods' however time saving is much more effi- 
cient via Eq. (p, since in fact this latter method corre- 
sponds to an average of the former over a large number of 
successful events. Both of them require additional mem- 
ory occupation for the storage of dimers and bonds that 
may change their state, but this does not seriously hinder 
simulations of large systems if compared with required 
computational time. 
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FIG. 4. Comparison of the temporal behavior of density of 
empty sites for an initial jammed lattice (solid line) and an 
initial empty lattice (long dashed line), for p = 0.5 The size 
of the system is L — 100. 



one dimensional lattices to form multilayer structures, 
that allow for both intra-layer and inter-layer diffusion, 
have been found to behave in a similar fashion.c3 



Our simulations show that also the RSAD of dimers in 
d = 2 displays a regime in which t -1 / 2 dominates, thus 
implying the recombination of one-site vacancies to be 
the main process, in analogy with the d — 1 case (we will 
turn to this point in Sec. |fV|) . This result has shown to be 
independent from the parameter values and initial con- 
ditions chosen for the simulation. Nevertheless in d = 2 
there is a stage of the covering process at which the dy- 
namics definitely changes. At that stage we observe that 
the annihilation rate slows down and correspondingly the 
filling slowly stops. The system is prevented to reach the 
complete coverage by the formation of isolated vacancies 
that are trapped at one site or along, closed loops. As 
firstly pointed out by Grigera et alx3 this seems to be 
a genuine property of the system, since independent of 
the lattice size. The results displayed in Fig. [l], where 
systems of different sizes approach the same asymptotic 
coverage 0(oo) ~ 0.998, support the idea of a dynami- 
cally jammed state related to local configurations_pf the 
system. We find, in agreement with the authors ho, that 
finite-size effects are reasonably negligible for L > 40. 



III. TIME COVERAGE 
A. Asymptotic coverage 

Mean field approximations are known to work very 
badly for RSAD of small particles in one dimension. 
Diffusional dynamics of dimers along the lattice can be 
described in terms of mobile vacancies, that annihilate 
through the deposition of a dimer when they meet on 
neighbor sites. As already mentioned, both exact solu- 
tions and numerical simulations show that the process of 
vacancies annihilation is dominated by random fluctua- 
tions of density rather than by average values, thus lead- 
ing asymptotically to a coverage described by a random 
walk process, 0(oo) — 6{t) ~ i -1 / 2 , opposed to the t~ x 
coverage estimated by mean field. Random walk process 
with regular exponent 0.5 concerns the dynamics of va- 
cancies, since particles suffer an anomalous,random walk 
with exponent < 0.5, as discussed in Ref.ta. 

Mean field is generally speaking expected to work bet- 
ter in higher dimensions. However studies on RSAD 
of different objects in d = 2 found that fluctuations 
of isahdfd vacancies still dominate relaxational dyjiarn- 
ics.Eja'Ej Investigations of Eisenberg and BaramEj on 
RSAD of monomers particles with nearest neighbors 
and next nearest neighbors exclusion have shown that 
9(oo) — 6{t) ~ £ -1 / 2 to the leading order in the asymp- 
totic regime. Systems made by superposition of several 



In Fig. g we show the coverage behavior for different 
choices of p. Different values of 9{oo) are reached, imply- 
ing that different choices of p and R can lead to different 
final coverages and can therefore be regarded as different 
quench rates of an initially annealed system. 
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FIG. 5. Temporal behavior of density of empty sites for 
infinite deposition rate (solid line). The long dashed line is 
a fit to the data according to Eq. (|E|). For comparison we 
also show the curve for p = 0.5 (dot-dashed line). The initial 
condition is 6* (0) = and the size of the system is L = 100. 
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FIG. 6. Temporal behavior of deferent kinds of vacancies 
for p = 0.2 (a), p = 0.5 (b) and p = 0.9 (c), for 0(0) = 0j: the 
solid line is the total density of empty sites, the dotted line is 
the density of one-site vacancies, the dashed line is the den- 
sity of sites formed by two-site vacancies and the dot-dashed 
line is the density of sites formed by non one-site vacancies. 
The size of the system is L — 100. 



B. System evolution 



Fig. |^ puts into evidence a somewhat counterintuitive 
feature, i.e. that faster coverages are obtained at lower 
deposition probabilities, showing that to this aim the pro- 
cess of vacancies generation is more important than depo- 
sition attempts themselves. In order to understand this 
behavior it is worth to consider the different processes 
that contribute to increase the system coverage. This 
can be done by describing the dynamics of the system on 
a wide time scale by a series of the form 



i=0 



(4) 



A fit of this curve to the results of numerical simulations 
shows that the terms beyond the fourth in the series can 
be neglected (cfr. Table ||) and is displayed in Fig. |[ As 
expected, the effect of an increased diffusion is to decorre- 
late the system, and thus to produce a faster coverage, as 
shown by higher weights for the i _1 term. The coverage 
depends on the initial condition only at early times. 
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FIG. 7. Temporal behavior of different kinds of vacancies 
for p = 0.2 (a), p = 0.5 (b) and p = 0.9 (c), for 0(0) = 0: the 
solid line is the total density of empty sites, the dotted line is 
the density of one-site vacancies, the dashed line is the den- 
sity of sites formed by two-site vacancies and the dot-dashed 
line is the density of sites formed by non one-site vacancies. 
The size of the system is L = 100. 

Fig. H displays the time coverage when 0(0) = 0. After 
an initial stage (1 < t < 10 2 ) in which deposition at- 
tempts dominate and the density progressively increases, 
the system approaches the same regime as for 9(0) = 9j. 

TABLE III. Coefficients of fit Eq. (|) for 0(0) = 0j. 



p 


X 


Xl 


X2 


2^3 


Xi 


0.2 
0.5 
0.9 


0.00131 
-0.00158 
0.00117 


0.093 
0.174 
0.249 


0.834 
0.602 
0.575 


2.932 
-2.916 
-4.801 


2.833 
3.016 
6.437 



TABLE IV. Coefficients of fit Eq. (g) for 0(0) = 0.; 



p 


2/o 


Vi 


V2 


J/3 


0.2 
0.5 
0.9 


0.00222 
0.00306 
-0.00010 


0.0509 
-0.0946 
-0.0033 


0.371 
0.925 
0.0896 


-0.446 
-1.484 
-0.1268 
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FIG. 8. Temporal behavior of empty sites density for 
p — 0.5 (a) and p — 0.1 (b), at different values of particle 
length k: k = 2,3,4,5 (from bottom to top in each graph). 
The initial condition is 8(0) = 8j and the size of the system 
is L = 100. 

A comparison of the two cases is made in Fig. ||| This 
time behavior can also be well described by curves of the 
form of Eq. (|]), obtaining the coefficients of Table |J. 
Coefficients relevant to the late stage process exhibit a 
similar trend with respect to those of Table |[ 

Another interesting situation, also often investigated 
in literature, is the Diffusion Limited Deposition (DLD). 
In this case p — 1 and R — oo, corresponding to the fill- 
ing of two-site vacancies with dimers as soon as they are 
created by diffusion. In this case only one-site vacancies 
exist, making the coverage process much more effective 
with respect to the case of finite R, as evidenced in Fig. ||. 

On the other hand the asymptotic coverage is higher 
for finite deposition rate, showing that the larger diffu- 
sion allowed becomes important at very high coverage 
values, since in that case the long time frozen config- 
urations are more easily removable. Thus, in analogy 
with Figs. [2] and |L different asymptotic coverages may 
be reached. Also for DLD a superposition of non ana- 
lytic powers of time fits the curve quite well in all the 
dynamical range: 



1 - 6{t) ~b + M~ 1/2 + ht- 1 + b a t 



-3/2 



(5) 



IV. VACANCY DYNAMICS 

As mentioned above, RSAD can be described in terms 
of reaction and diffusion of vacancies on the lattice. In 
the general case of k— mers with diffusion coefficient D 
and deposition rate R, one deposition corresponds to the 
annihilation of k vacancies on contiguous sites. 

If one adopts this point of view the system behavior ap- 
pears to be ruled by the dynamics of vacancies. In order 
to give a full account of this dynamics, knowledge of va- 
cancies correlation functions is necessary. Nevertheless a 
qualitative understanding of some of the main dynamical 
features can be obtained also through the investigation 
of the sole vacancy population. 



A. Initial jammed state 

Fig. U compares the coverage behavior with the change 
in time of vacancy population for different values of p. 
Solid and dotted lines represent respectively the total 
fraction of empty sites, 1 — 0, and of one-site vacancies, 
pi . At early and intermediate times these two quantities 
display some differences, more visible at lower values of 
p where formation of many-site vacancies is enhanced by 
diffusion. For a jammed initial state, an increase in two- 
site vacancies takes place in the early stage of the process 
and, for low values of p, can last for several decades both 
in density (^2) and time. In fact, two dimensional ge- 
ometry allows for the formation of many different kinds 
of voids; besides vertical and horizontal two-site vacan- 
cies, there can be L many-site vacancies with different 
orientations, square vacancies and so on. 

TABLE V. Jamming densities of k— mers. 



0.8470 
0.8103 
0.7877 



TABLE VI. Coefficients of fit Eq. (§) for fc-mers when 
?(0) = 9 j. 



where b„ ~ 0.0016, h ~ 0.060, b 2 
—0.137 as shown in Fig. 
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0.399 
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10.89 


5 (p = 0.5) 


0.00228 
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-1.774 


0.8537 
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2 (p = 0.1) 


0.00126 
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37.44 


5 (p = 0.1) 
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FIG. 9. Temporal behavior of density of one-site vacan- 
cies (solid lines) and of sites formed by non one-site vacancies 
(dotted lines) for p = 0.1 and k = 3 (a), k — 4 (b), k = 5 (c). 
The initial condition is 6(0) = 8j and the size of the system 
is L = 100. 

A comparison of P2 with the overall non-one-site va- 
cancy density, 1 — p± — 8, is shown in Fig. ||. We observe 
that the abundance of the latter does not vanish for large 
times, reflecting the configurational richness of the two- 
dimensional process with respect to the unidimensional 
one where p2 and 1 — p\ — 8 coincide asymptotically. 

If one expands p2 and 1 — p\ — 8 in powers of i -1 / 2 
one finds, as for the coverage, that larger contributions 
from mean field terms appear at lower p values in the 
asymptotic expansions: 



P2(t) 



i=0 



■i/2 



■i/2 



(6) 



(7) 



i=0 



The coefficients of the best fit of Eqs. ^ and (0) to 
the numerically determined curves for p2 and 1 — p \ — 8 
are reported respectively in Table III and Table |Vl 



B. Initial empty lattice 

An initially empty lattice produces early and inter- 
mediate stages in vacancy dynamics different from the 
jammed case. Whereas the late stage is still dominated 
by recombination of one-site vacancies, an initial stage 
exists in which the abundance of this kind of defects in- 
creases monotonically up to reach about the same density 
of many-site vacancies, as displayed in Fig. 0. On the 
other hand, right at that stage the latter's population 
starts to be less than that of one-site vacancies. Inter- 
estingly single empty sites rule the dynamics for t > 10 2 , 
that is, as already explained, the same time range in 
which density no longer depends on initial conditions 
(cfr. Fig. ||). It can be noted that two-site vacancies 
grow at pace with one-site ones up to a maximum, where 
they start to decay at the same rate of larger ones. The 
location of this maximum as well as the intersection be- 
tween 2p2 and 1 — pi — 8 tend to move to early times 
when p increases, suggesting that enhanced deposition 
rate destroys many-site vacancies more efficiently. The 
functional forms that fit well p2 and 1 — p± — 8 are the 
same as those for the initial jammed lattice case (cfr. 
Eqs. (§) and §)). 



V. RSAD OF K-MERS 

We now address the problem of RSAD of k— mers on a 
square lattice, which is a generalization of the model de- 
scribed in Sec. || that corresponds to k = 2. Since for in- 
creasing k the simulations of the dynamics becomes more 
and more time consuming and finite-size effects more 
important, we concentrate on low values of k, namely 
ft = 3, 4, 5. 

We have simulated the RSAD dynamics always start- 
ing with an initial jammed state and then switching on 
the diffusive part of the process. We checked that the 
jamming densities obtained as initial conditions after the 
pure RSA process were in aweement with the values 
found earlier by Bonnier et alE3 It turns out that there 
were no discrepancies with their results within the nu- 
merical precision. The jamming densities obtained are 
summarized in the Table |v|. 

In Fig. H we compare the temporal evolution of the 
ft— mers density for two values of p. 



TABLE VII. 
tion. 



Crossover time for one site vacancy domina- 



te 



80 
400 
1680 
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The major aspect to underline is the very sluggish dy- 
namical behavior for increasing k, which is a common 
feature with the one dimensional case: the probability 
for k vacancies to meet and annihilate decreases when k 
gets larger. It can be noted that, as for dimers, a dynam- 
ical saturation at long times exists, although our data are 
less clear for k = 5, since the time at which this plateau 
sets in is an increasing function of k. Such a saturation is 
absent in d = 1, where the density follows apSingle power 
law for long times also for larger values of fc.Eil It appears 
therefore to be a feature typical of the two dimensional 
geometry and independent of the particle size, related to 
the complexity of the configurational space accessible to 
the system. A proof of that is that the curves in Fig. || 
can be fitted by the same functional form adopted for 
dimers (Eq. (Q)) with the parameters of Table |vi| It is 
seen that the i -1 / 2 term is more important at larger k, 
in agreement with a slow down in the coverage process. 

Arguments based on diffusive reactionscll have shown 
in the past that a critical dimension d c — 2/(k— 1) should 
exist above which fluctuations would be negligible and 
mean field behavior (1 — 6 ~ t -1 '^ 1 ') would dominate. 
This would imply that for d = 2 the mean field law should 
work for k > 3; nevertheless we have found no evidence 
of that, at least for k — 3. In order to draw such a con- 
clusion for larger k, longer runs would be required, since 
a change in the dynamical behavior could occur at very 
large times. 

Finally we have compared the data for the density 
of one-site vacancies and of many-site vacancies, for 
k = 3,4,5 and p = 0.1 (Fig. |J). At variance with the 
case k — 2, many-site vacancies tend now to prevail for 
a long time and disappear more slowly for larger k; the 
time behavior of their density can be well approximated 
by a single power law: 



particular, the existence of dynamically frozen configura- 
tions, that only depend on the local environment, slows 
down the asymptotic dynamics and does not allow the 
system to reach the complete coverage. Different choices 
of deposition and diffusion rates and probabilities can 
lead to different final values of the asymptotic coverage. 
Simulation shows that the behavior of coverage strongly 
deviates from the mean field prediction. Deviations from 
mean field are less pronounced when the diffusion proba- 
bility is higher. Furthermore 9 is observed to depend on 
the initial conditions only in the early stage of the dy- 
namics. Starting from the jammed state correlations are 
higher and the dynamics evolves more slowly. Vacancy 
dynamics is dominated not only by single-site vacancies, 
as in d — 1, but also by larger vacancies. 

It was also shown that the dynamical behavior is 
severely slowed down with the increase of the length k of 
the particles. We do not have evidences of a mean field 
behavior, in spite of the fact that it would be expected 
to hold for k > 3. We have also pointed out that the di- 
mensionality plays an important role in reaction-diffusion 
models with the result of a rich dynamical behavior in the 
investigated quantities. 

As an open possibility for the future, we think that in 
order to better understand the dynamics of these models 
it would be useful to vary the kind of relaxation mecha- 
nisms, investigating the possibility of occurrence of col- 
lective effects in d = 2 and d = 3. 
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l-pi(i)-0(t)~f 



-fa 



(8) 



with /3 3 = 0.94, /? 4 = 0.75 and (3 5 = 0.61. Moreover, the 
crossover time t c after which one-s ite v acancies dominate 
increases with k, as shown in table VII of t c values, where 
we note that for higher k values diffusive regime sets in 
at later time. 



VI. SUMMARY AND CONCLUSIONS 

We studied irreversible deposition and diffusion of 
dimers and k-meis on a square lattice, using a driven al- 
gorithm which allowed to decrease the CPU time needed 
for Monte Carlo simulations. Both time coverage behav- 
ior and vacancy dynamics have been analyzed in detail 
as a function of initial conditions and deposition rate and 
probability. A wider variety of dynamical regimes with 
respect to the one dimensional case has been observed, 
reflecting the complexity of the configurational space ac- 
cessible to dimers in d = 2 with respect to d = 1. In 
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